Neutrino Flux from Cosmic Ray Accelerators in the Cygnus Spiral Arm of the Galaxy 
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Intriguing evidence has been accumulating for the production of cosmic rays in the Cygnus region of 
the Galactic plane. We here show that the IceCube experiment can produce incontrovertible evidence 
for cosmic ray acceleration by observing the neutrinos from the decay of charged pions accompanying 
the TeV photon flux observed in the HEGRA, Whipple, Tibet and Milagro experiments. Our 
assumption is that the TeV photons observed are the decay products of neutral pions produced by 
cosmic ray accelerators in the nearby spiral arm of the galaxy. Because of the proximity of the 
sources, IceCube will obtain evidence at the 5a level in 15 years of observation. 



Evidence may be emerging for a cosmic accelerator in 
the Cygnus spiral arm. The observation of the Cygnus 
region by the HEGRA IACT-system has allowed the 
serendipitous discovery of a TeV 7-ray source [l[, with 
an average flux ~ 3% of the Crab Nebula [2j . The analy- 
sis of the total 278.3 hours of observations performed in 
two periods from 1999 to 2002 (120.5 hours from 1999 to 
2001 [H and 157.8 hours during 2002 [H) has revealed the 
presence of a steady (and possibly extended) TeV source, 
with hard injection spectrum. 

The excess significance of the TeV source is 7.1c and it 
appears extended at more than 4cr level with a morphol- 
ogy which is suitably described by a Gaussian profile. 
The source is termed TeV J2032+4130 after the position 
of the center of gravity and its extension (Gaussian la 
radius) is 6.2'(±1.2g tat ± 0.9g ys ). Especially intriguing is 
the possible association of TeV J2032+4130 with Cygnus 
OB2, a cluster of more than 2700 (identified) young, hot 
stars with a total mass of ~ 10 4 solar masses Q. At a 
relatively small distance (« 5000 light years) to Earth || , 
this is the largest massive Galactic stellar association. 

The observed hot spot has no clear counterpart and the 
spectrum is not easily accommodated with synchrotron 
radiation by electrons. The difficulty to accommodate 
the spectrum by conventional electromagnetic mecha- 
nisms has been exacerbated by the failure of CHANDRA 
and VLA to detect X-rays or radiowaves signaling accel- 
eration of any electrons [6j . The two most plausible mod- 
els to explain the 7-ray signal are: (a) A proton beam 
(accelerated in the stellar winds of Cygnus OB2 [J, @], 
or else in the wind nebulae of an undetected nearby pul- 
sar Q) interacting with a molecular cloud to produce 
pions that are the source of the gamma rays. Proton 
acceleration to explain the TeV photon signal requires 
only 1% efficiency for the conversion of the energy in the 
stellar wind into cosmic ray acceleration, (b) The TeV 
gamma rays can also originate in the photo-deexcitation 
of ultra-relativistic nuclei (Lorentz factor ss 10 6 ) that are 
themselves the photo-disintegration products of heavier 
nuclei broken-up in the bath of intense UV photons from 
the Lyman a emissions of hot stars Q . As in the proton 
beam model, the required power density for acceleration 



of nuclei is 2 orders of magnitude smaller than the kinetic 
energy budget of the entire association. 

An additional set of observations performed during 
1989-1990 by the Whipple Observatory [§] has been re- 
cently reanalyzed in the light of the HEGRA data. These 
confirm an excess in the same direction as J2032+4130, 
although with considerably larger average flux (~ 12% of 
the Crab), above a peak energy response of 0.6 TeV. The 
statistical significance of the signal is only 10% smaller 
with selection of events above 1.2 TeV. However, the 
large differences between the flux levels cannot be ex- 
plained as errors in estimation of the sensitivity of the 
experiments since they have been calibrated by the si- 
multaneous observations of other TeV sources. More re- 
cently, data taken with the Whipple Observatory dur- 
ing 2003-2005 have been reported The analysis of 
the latest dataset reveals a TeV hot spot (integral flux 
~ 8% of the Crab) that is displaced about 9 arcmin- 
utes to the northeast of the TeV J2032+4130 position. 
A re-analysis [ll[ of a 10-hour VLA mosaic exposure to- 
wards TeV J2032+4130 (based on the alternative source 
hypothesis) allowed the detection of a weak, predom- 
inantly non-thermal, shell-like supernova remnant-type 
object (with location and morphology very similar to the 
HEGRA source) that can be the cosmic ray engine pow- 
ering the OB association. 

Very recently, the Milagro Collaboration reported an 
excess of events from the Cygnus region at the 10.9<r 
level [l2j. The observed flux within a 3° x 3° window 
centered at the HEGRA source is 70% of the Crab at 
the median detected energy of 12 TeV, and has a dif- 
ferential spectrum oc E~ 2e . Such a flux largely exceeds 
the one reported by the HEGRA Collaboration, implying 
that there could be a population of unresolved TeV 7-ray 
sources within the Cygnus OB2 association. 

The Milagro Collaboration also reported a new hot 
spot, christened MGRO J2019+37, at right assension 
= 304.83° ± 0.14 stat ± 0.3 sys and declination = 36.83° ± 
0-08 s tat ± 0.25 sys [l2j. This new unidentified source is 
observed with statistical significance > 6a above the av- 
erage diffuse 7-ray emission in the region. A fit to a cir- 
cular 2-dimensional Gaussian yields a width of 0.32±0.12 
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degrees, which for a distance of 1.7 kpc suggests a source 
radius of about 9 pc. For a differential spectrum oc E~ 2S 7 
the brightest hotspot in the Milagro map of the Cygnus 
region represents a flux of 1 Crab above 12.5 TeV. In- 
terestingly, the Tibet AS-gamma Collaboration has ob- 
served a cosmic ray anisotropy from the direction of 
Cygnus, which is consistent with Milagro's measure- 
ments [13] . Unfortunately, the Tibet array has very little 
power to distinguish how much of the anisotropy should 
be attributed to 7-rays and how much, if any, to baryons. 

The brightest Milagro hot spot is located outside the 
OB association. However, the 7-ray signal is found to 
trace the gas density distribution in the region. The 
model proposed is that of a cosmic ray beam, perhaps 
powered by a millisecond pulsar, which interacts with a 
molecular cloud positioned a few degrees to the south- 
east of the OB star cluster. If the 7-ray emission from 
MGRO J2019+37 originates in ir° decay, it is necessarily 
accompanied by a flux of high energy neutrinos emerging 
from the tt population. In this paper we discuss in de- 
tail the prospects to observe such a flux with the IceCube 
neutrino telescope. 

The pion spectrum resulting from collisions of the 
ultrarelativistic protons on the molecular cloud is ex- 
pected to obey a modified Feynman scaling in the cen- 
tral rapidity region, dN w /dE^\E p ~ C (E p ) / E^ , where C 
may be growing as some power of In E p [l4j . For given 



dE scale with these fractions, and we arrive at [if 



En < 0.08 E, 
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we may convolve with a proton spec- 



trum typical of Fermi engines, dN p /dE p cx E r , 
tain the pion spectrum ( 151 ] 
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where C{E 7r ) is generically a function which grows as 
a power of ln-E^, falling to zero at the cutoff E n = 
0.08 E p<inax . Since 7r°'s, 7r + 's, and 7r _ 's are made in equal 
numbers, one expects two photons, two v^s, and four 
v^s per 7T°. Gamma rays, produced via tt° decay carry 
one-half of the energy of the pion. Each ir~ decays to 
3 neutrinos and an electron, n~ — ► pTV^ — > Vfjj^v^ . 
The electron radiatively cools through interactions with 
the gas and the ambient magnetic and radiation fields. 
Typically e~ synchrotron emission extends from radio 
frequencies to X-rays. The average neutrino energy from 
the direct pion decay is {E v ) v = (1 — r) E w /2 ~ 0.22 E n 
and that of the muon is (£' M > 7r = (1 + r) E^/2 ~ 0.78 E„, 
where r is the ratio of muon to the pion mass squared. 
Now, taking the from muon decay to have 1 /3 the en- 
ergy of the muon, the average energy of the from muon 
decay is (E^^)^ = (1 + r)E K /6 = 0.26 E^. Similar con- 
siderations apply for the charged conjugate process. For 
simplicity, hereafter we consider that all neutrinos carry 
one quarter of the energy of the pion. The energy-bins 
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where tt denotes any one of the three pion charge-states. 

Whereas the details are complex and predictions can 
be treacherous, it is clear that the astrophysical ambi- 
guities far outweigh the details associated with the par- 
ticle physics, and hence it is safe to assume that iden- 
tical fluxes of 7-rays and are produced. Terrestrial 
experiments have shown that i/„ and v r are maximally 
mixed with a mass-squared difference ~ 10 _3 eV 2 , and 
that |(i/ e |i/ 3 )| 2 is nearly zero [ijj ■ Here 1/3 ~ {Vfj, + v T )/ y/2 
is the third neutrino eigenstate. This implies that any 
initial flavor ratio having tu e = 1/3 will arrive at Earth 
with ratios u) e : u)^ : u> T = 1 : 1 : 1. Thus, there is a fairly 
robust prediction that the initial flavor ratios of 1 : 2 : 
given in Eq. ([3]) would arrive at Earth democratically dis- 
tributed, i.e., 1 : 1 : 1. From these remarks, one finds a 
nearly identical flux, 
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for each of the three neutrino flavors a = e, fi, r [l8|. 

The Antarctic Muon And Neutrino Detector Array 
(AMANDA) using natural 1 mile deep Antarctic 

ice as a Cerenkov detector, has operated for more than 5 
years in its final configuration: 19 strings instrumented 
with 680 optical modules. IceCube [2(|, the successor ex- 
periment to AMANDA, is now under construction. It will 
consist of 80 kilometer- length strings, each instrumented 
with 60 digital optical modules (DOM) spaced by 17 m. 
The deepest module is 2.4 km below the surface. The 
strings are arranged at the apexes of equilateral trian- 
gles 125 m on a side. The instrumented (not effective!) 
detector volume is a full cubic kilometer. A surface air 
shower detector, IceTop, consisting of 160 Cerenkov de- 
tectors deployed over 1 km 2 above IceCube, augments the 
deep-ice component by providing a tool for calibration, 
background rejection and air-shower phy sics. The angu- 
lar resolution for muon tracks ~ 0.7° [2l| allows a search 
window of 1° x 1°. Construction of the detector started 
in the Austral summer of 2004/2005 and will continue 
for 6 years, possibly less. At the time of writing, data 
collection by the first 22 strings and 52 IceTop stations 
has begun. 

The event rate from a source located at given declina- 
tion can be calculated from the knowledge of the so-called 
neutrino effective area. This parameter strongly depends 
on the energy due to the almost linear increase of the 
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FIG. 1: The IceCube effective area for muon track recon- 
struction as a function of the neutrino energy. Two curves 
are shown indicating the trigger level and the recent estimate 
with very conservative quality cuts. For comparison the ef- 
fective area of AMANDA II is also shown. 



FIG. 2: Differential energy distribution of event rates from 
MGRO J2019+37. The different curves indicate rates com- 
puted using the effective area (a) of the trigger level condition, 
(6) of the conservative quality cuts, (c) of the AMANDA II 
array. 



neutrino interaction cross-section and to the logaritmic 
rise of the muon range at very high energies. Moreover 
it accounts for the probability of absorption of neutrinos 
during their propagation through the Earth, where they 
may disappear due to a charged current interaction. It 
should also account for the detection and reconstruction 
efficiency. Preliminary studies by the IceCube Collabo- 
ration [23] , using data collected with the 9 In-Ice strings 
and the 16 IceTop stations, show a good agreement with 
the detector performance [20| . The IceCube simulation 
software is currently under active development and ex- 
trapolations of the 9-string results are undeway. In par- 
ticular, a conservative estimate of the effective area for 
reconstruction of muons tracks, with severe cuts to ob- 
tain a good angular resolution and rejection of misrc- 
constructed muons has been estimated in [22]. Such an 
effective area is shown in Fig. [Tj together with the trig- 
ger level effective area (defined as 8 DOM threshold) of 
IceCube and the effective area of the AMANDA final 
configuration. 

Equipped with the effective area shown in Fig. [TJ it is 
straightforward to calculate lower and upper limits on 
the muon neutrino event rate from the Cygnus region 
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Using the neutrino flux given in Eq. and convolving 
it with the IceCube neutrino effective area we foresee 
an event rate of muon neutrinos with E v > 1 TeV of 
1.1 yr^ 1 < dNs/dt < 4.1 yr _1 . The energy distribution 
of such events is shown in Fig. [2] 

We now turn to the estimate of the background. For 
the atmospheric neutrino flux, arising from the decay of 
pions and kaons produced in cosmic ray interactions with 
the air molecules, we adopt the estimates of Ref. [23|. We 
obtain the number of expected muon tracks from atmo- 
spheric neutrinos as in Eq. JTJ) , using the atmospheric 



neutrino flux integrated over a solid angle of 1° x 1° width 
around the direction of the MGRO J2019+37 (zenith an- 
gle 9 = 53.2°). We obtain an expected background of 
atmospheric tracks, 1.2 yr _1 < dNs/dt < 5.5 yr" 1 . 

These event rates are based on a conservative estimate 
of the level of detail to which we currently understand the 
detector performance. However, since this understanding 
will improve over time one expects the systematic errors 
to decrease to the levels projected in the baseline design 
reported in [2fj| . Therefore, to determine the discovery 
reach we employ the semianalytical calculation presented 
in [24| based on a full Monte Carlo simulation using these 
projected baseline detector properties, with quality cuts 
referred as level 2 cuts [2(| . For a muon energy threshold 
of 100 GeV and minimum track length of 300 m, the 
expected rate of induced tracks is dNs/dt ~ 3 yr _1 , 
with a background of dNs/dt ~ 2.5 yr" 1 . Hence, after 
15 yr of operation, the (total) detection significance, 
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is expected to be at discovery level. 

We now verify that our results are consistent with ex- 
isting data. Very recently, the AMANDA II data col- 
lected during 2000 - 2004 (with a lifetime of 1001 days) 
was analyzed to set new limits on the neutrino fluxes from 
point sources (circular bin size varying between 2.25° and 
3.75° depending on declination) [25j. For the effective 
area shown in Fig. [TJ the expected background from at- 
mospheric neutrinos {E v > 100 GeV) in a 1° x 1° window 
is dN^/dt ~ 0.22 yr- 1 . The AMANDA Collaboration 
reported an expected atmospheric neutrino background 
from the direction of J2032+4130 of 6.8 events, in com- 
plete agreement with our calculations. The experiment 
observed 7 events pointing towards the same direction of 
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the sky, leading to a 90% CL upper limit, 
dF 

El—^- = 1.1 x 1CT 10 TeVcm" 2 s" 1 . (6) 

Thus, the sensitivity reach at AMANDA cannot probe 
the predicted flux given in Eq. ©. 

In summary, by observing the neutrinos from the de- 
cay of charged pions accompanying the recently detected 
7-rays with the Milagro experiment, the IceCube neu- 
trino telescope will produce incontrovertible evidence for 
cosmic ray acceleration in the Cygnus spiral arm. In this 
paper, we have discussed in detail the sensitivity reach 
of IceCube to the brightest hot spot, MGRO J2019+37. 
Contributions from TeV 7-ray hot spots (still hidden) 
within the Cygnus OB2 association will certainly en- 
hance the signal [27|. Moreover, as suggested recently 
by the Milagro Collaboration [Tj|], to smoothly match 
EGRET data in the 100 MeV energy region, the spectra 
from all the sources in the Cygnus region should have 
a break (perhaps because of absorption effects) and be 
harder than E~ 2 6 at lower energies. As an illustra- 



tion, we have estimated the expected event rate from 
MGRO J2019+37, using the semianalytical calculation 
given in [24[ and assuming a spectrum cx E~ 2A , 

dNs/dt ~ 9.0 yr" 1 . (7) 

For a background of 2.5 yr _1 events, this implies that 
IceCube will attain a 5er discovery reach in 2 years of 
operation! 
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